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Abstract 

Chiral quark models offer a practical and simple tool to describe covariantly 
both low and high energy phenomenology in combination with QCD evolution. 
This can be done in full harmony with chiral symmetry and electromagnetic 
gauge invariance. We review the recently proposed spectral quark model where 
all these constraints are implemented. 

1. INTRODUCTION 

Light cone hadron properties provide a non-trivial playground to test ideas and methods on genuine non- 
perturbative QCD phenomenology. As we have heard from S. Brodsky in this workshop this information 
is conveniently encoded in LC wave functions Q. They are measured through high energy inclusive, 
such as Deep Inelastic Scattering (DIS), or exclusive processes, which with the help of QCD renormal- 
ization group evolution can be related to low-energy matrix elements. In the particularly interesting case 
of the pion one expects chiral symmetry to impose certain constraints. The specific and general way 
how this is done seems unclear at the moment, because by asking about the quark content of a hadron 
we are actually in a situation where both hadrons and quarks coexist as degrees of freedom, so standard 
effective field theory tools like Chiral Perturbation Theory do not apply directly beyond predicting 
finite pion mass corrections. On the other hand, fundamental lattice approaches naturally formulated in 
Euclidean space can only obtain low moments of structure functions (see e.g. the talk of P. E. L. Rakow 
and Ref. [3]) and distribution amplitudes, which are not testable directly through experiment. As B. 
van de Sande presented in his talk, QCD transverse lattice approaches can be directly formulated on the 
light-cone and although promising results for the pion and other mesons exist already [4| (for a review 
see e.g. Ref. J5 ]) they do not easily accommodate for chiral symmetry. The Schwinger-Dyson approach 
in its Euclidean formulation and its implications [6| has been reviewed by C. D. Roberts. 

In order to gather some theoretical insight on the role of chiral symmetry in high energy inclu- 
sive Q or exclusive (8 1 processes we remain at the more modest level of the so-called quiral quark 
models reviewed in Ref. [9]. This generic name stands for phenomenological relativistic quantum field 
theories where chiral symmetry is spontaneously broken incorporating Goldstone's theorem. One of 
the advantages of such an approach is that it is covariant, so we avoid from the beginning the difficult 
marriage between light-cone quantization and chiral symmetry. In this talk we review salient features of 
recent work done by us on a new chiral quark model: the spectral quark model (SQM) [ 103^D[^). This 
model shares many good aspects of other chiral models, particularly NJL type models, but improves oth- 
ers at a very low computational cost and, up to now, brings very encouraging phenomenological success. 

QCD is a theory of quarks and gluons as elementary degrees of freedom, so both contribute to 
the total momentum carried by a hadron. The separation between gluons and valence and sea quarks, 
however, is both scale and renormalization scheme dependent which can be worked out explicitly within 
perturbative QCD. On the other hand, chiral quark models which are supposedly a low energy approx- 
imation to QCD, obviously saturate the momentum sum rule, since they contain no explicit gluonic 
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degrees of freedom. This poses the natural question: is there any scale in QCD at which the momentum 
fraction carried by the (valence) quarks is exactly one ? 

According to the phenomenological QCD analysis undertaken by the Durham group over a decade 
ago HI 31 . the momentum fraction carried by the valence quarks is 0.47 % at the scale /j 2 = 4GeV 2 , e.g., 

for 7T + , 

V 1 (fj,) = (x(un-u n + d 7r -d 7r )) = 0A7±0.02 at i? = 4GeV 2 . (1) 

where u^,u^ , d^ and d n are parton distribution functions (PDF) in the pion. Assuming LO perturbative 
QCD evolution one gets 
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where 7f s / 2 A) = 32/81 for N F = N c = 3. We take for concreteness A QC d = 226 MeV, which for 
fi = 2GeV yields a = 0.32. Downward LO evolution yields that for the reference scale, (j,q, 

y 1 (/i ) = l, M = 313i 2 °MeV (3) 

Although this seems a rather low scale, one may still hope that the expansion parameter a(fio)/2ir ~ 
0.34±0.04 makes perturbation theory meaningful. A NLO analysis confirms this, at first glance "illegal", 
expectation. This is the natural scale where all observables are defined in the quark model. If we want 
to compute observables at higher scales fj, > [Iq one may use QCD (perturbative) evolution using the 
quark model as an initial condition. This way we generate the missing (perturbative) gluon components 
explicitly in the form of QCD radiative corrections. 



2. THE SPECTRAL QUARK MODEL 

The subject of regularization in connection with high energy hadronic properties in chiral quark models 
has always been tricky and frustrating in the past, particularly when chiral and electromagnetic invari- 
ance, anomalies and factorization were demanded simultaneously (See Ref . [ 9 1 and references therein 
for a comprehensive discussion). The SQM relies on a spectral regularization of the chiral quark model 
proposed in Ref. [ 10 1 and extensively developed in Ref. [ 11 1 and is based on a generalized Lehmann 
representation for the quark propagator 1 

S{p) = / duj- = / dw 2 2 = 1 ^77~^v W 

Jc p-uj Jc P - uj z p- M{j> 1 ) 

where C is a general contour to be determined. For obvious reasons uj is called the spectral quark mass. 
Chiral and flavour gauge invariance on the relevant vertex functions are obtained (up to transverse terms) 
by means of the gauge technique of Delbourgo and West [14] which is nothing but minimal coupling 
at the spectral quark level. According to the analysis of Refs. fl"0l IXTIl . the proper normalization and 
finiteness of hadronic observables are achieved by requesting a set of spectral conditions for the moments 
of the quark spectral function, p(u), namely 

p n = J duu n p(u) = 5 n0 , for n = 0, 1, 2, 3, ... (5) 

Physical observables turn out to be proportional to the inverse moments, 



= J dujuj k p(uj), for k= 1,2, 3, 



(6) 



'Note that an approach similar in spirit was described long ago by Efimov and Ivanov 1151 . 
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Fig. 1 : The dependence of the quark mass M (top) and the wave function normalization Z (bottom) on the Euclidean momentum 
Q obtained from the meson dominance model with ds = 5/2. Lattice data from Ref. 1161 . 



as well as to the so called "log moments" , 

p' n = duo \og(u> 2 / p 2 )uj n p{uj) = I du log(u> 2 )co n p(uj) , for n = 2, 3, 4, ... 



Note that the conditions remove the dependence on the scale p in @, thus we can drop it from the 
log. Using these conditions one may prove a number of important features, like fulfillment of chiral 
anomalies and factorization of form factors, i.e. power like behaviour, in the high momentum limit. The 
main assumption in Eq. © is that of analyticity in the complex plane but not of positivity, since then 
even moments would not vanish for a positive spectral function. Analyticity enables direct calculations 
of hadron properties in the Minkowski space, an extremely convenient circumstance when doing light- 
cone physics. The analyticity assumption is also necessary in models formulated in Euclidean space, 
although the practical continuation to Minkowski may become numerically messy. 

Although many results of our model do not depend on a particular choice of functions fulfilling 
the spectral conditions ©, a particular realization embodies more predictions. The standard approach to 
quark models is to assume a "reasonable" form for the quark propagator and work out hadronic properties 
from there. This poses the natural and difficult problem as to what is a reasonable spectral function. 
Out of ignorance we take a radically different approach, namely we deduce the spectral function from 
phenomenology. Actually, py(u>), can be uniquely deduced from the pion electromagnetic form factor, 
which for simplicity is assumed to be of the vector dominance monopole form. The net result ifTTl is 
the meson dominance (MD) model explicit realization which fulfills the above spectral conditions and is 
given by 

( \ 1 1 1 ( , 1 16(ds-l)(rf s -2y 3 

PV[UJ> 2iTiuj(l-Auj 2 /M v ) d v' PS[UJ> 2ni M|(l - 4u 2 /M^s ' {,) 

The preferred values are Ay = ds = 5/2 and My = 776MeV corresponds to the p-meson mass. The 



contour C encircles the branch cuts, i.e., starts at — oo + iO, goes around the branch point at —m p /2, 
and returns to — oo — iO, with the other section obtained by a reflexion with respect to the origin 1111 . 
Straightforward calculation yields 



M(p 2 ) U v p 2 \M*-4p* 



2 \ d S 



M(0) My ( Ml 



M£-4p 2 



(8) 



1 



Remarkably, when M(p 2 ) = p 2 then Z(p 2 ) = so that the quark propagator has no pole. Instead it has 
cuts at p 2 = My j 4. These feature is usually called analytic confinement and the rationale for it has been 
described in detail in Ref. [11 1. A fit to the recent QCD lattice simulation in Landau gauge [ 16] yields 
M (0) = 303 ± 24 MeV and M s = 970 ± 21 MeV, with the optimum value of x 2 per degree of freedom 
equal to 0.72, yielding the impressive agreement of M(p 2 ) shown in the Fig. Q. As we see Z(p 2 ) is not 
as good, so there is obviously still room for improvement. 

One of the main appeals of the SQM is its computational simplicity. Using the above spectral 
conditions one can easily compute many properties by doing standard one loop calculations. For instance, 
the quark condensate for one flavour yields 2 . 

/- \ N c f, , f 2, 3 ( \MD d V MlN c M(0) 

{qq) = "4^ J dul ° g{uJ )U P{U) = nds-l)(ds-2)M 2 n 2 - (9) 

With the MD spectral density, ps, given by Eq. @ this equation becomes an identity. The corresponding 
value of the quark condensate is (qq) = — (243. Olg § MeV) 3 at the model scale po, which after LO 
evolution yields (qq) = — (301.0±?;J MeV) 3 at p = lGeV 0. The vacuum energy density is scale 
independent can be also be computed to give 

5 = -^/^log(^W-) M = D -^ = —(217 MeV) 4 . (10) 

for three flavors, Nf = 3. These values are reasonable according to current QCD sum rules estimates 
= -(243.0 MeV) 3 (at p = lGeV) and B = -(224l^MeV) 4 (at any scale) HD. 



3. PION PROPERTIES 

The pion weak decay constant can be worked out to give 



fl = -^J ^log( W V/^) M =^ = (87MeV) 2 (11) 

for My = 776MeV in excellent agreement with value of f n = 88MeV in the chiral limit obtained from 
ChPT \ 2\. The pion charge form factor is given by 



F v (t) = Fr(t) = J du>p(u)u 2 jf dxlog[co 2 + x(l-x)t]^ Q2+ v M2 (12) 

with the built in vector meson dominance explicitly displayed. This form describes the data remarkably 
well with My ~ 730MeV up to Q 2 ~ 1.6GeV 2 IT9I . 

A calculation along the lines of Refs. [ 8 , 9 1 produces the following identity between the light-cone 
pion wave function and the unintegrated pion distribution function 

U \ ( u \ Nc I a ( \ ^ MD 3M v 



2 The upperscript MD means evaluated in the Meson Dominance Model, Eq. jTJ. 




Fig. 2: The pion distribution amplitude (PDA) in the chiral limit evolved to the scale Q 2 — (2.4GeV) 2 . The two values for the 
evolution ratio r = a{Q) / a(Qo) reflect the uncertainties based on an analysis 1201 of the CLEO data 1181 . We also show the 
unevolved PDA, (^(2;, Qo) = 1, and the asymptotic PDA, ip 1T {x,oo) — 6x(l — x). 



where we remind that the identity holds in our model at the low-energy scale, /irj, of the model. In 
the MD the average transverse momentum squared is equal to (kj_) = My/2 = (544 MeV) 2 . After 
integration over d 2 k± the result for the PDA and PDF holds, 

q(x) = if(x) = 1 for fj, = fiQ. (14) 

The trend of both PDA and PDF being one has also been observed in transverse lattice calculations |4). 
The Gegenbauer evolution taking into account the uncertainties discussed by A. Bakulev in his talk [20| 
produce the PDA depicted in Fig. ©. In addition one may deduce after DGLAP evolution of the PDF 
and Gegenbauer evolution of the PDA the following relation 

6x(l — x) 

where the explicit analytic expression for the scale independent kernel, K(x,y) is given in Ref. [8]. This 
is a remarkable equation since it relates an inclusive process to an exclusive one. Using this equation 
one can regard the PDA measured at CLEO as a prediction in terms of the PDF parameterizations of 
Durham Ref. [13|. The QCD DGLAP-evolved PDF provide an impressive description of the Durham 
parameterization 1131 . The comparison to Drell-Yan data at Q 2 = 16GeV 2 as well as the extension to 
non-skewed generalized partem distributions can be found in the contribution by WB in this workshop 
(see also Ref. [ 12]). Using Eq. (fT3t one gets the following estimate for the leading twist contribution to 
the pion form factor at LO 

Q 2 J WQ) _ f 1 dx ^Q) = L25 ± ai0 (16) 



1= f 1 dyK(x,y)V^{y,n) (15) 
Jo 



2f n Twist-2 J 6x(1 - x) 

The experimental value obtained in CLEO [ 18] for the full form factor is Q 2 F 7 * i7Ty (Q) /2f w = 0.83 ± 
0.12 at Q 2 = 5.8GeV 2 . Considering that we have not included neither NLO effects nor an estimate 
of higher twist contributions, the 2-sigma discrepancy is not unexpected. At k± = the following 
conditions are satisfied, 

= =?(M)- (17) 

The first equality corresponds to the proper normalization of the anomalous decay 7r° — > 27. Other 
anomalous vertices are preserved as well ATI . 



To summarize, the SQM in conjunction with QCD evolution provides an appealing and compu- 
tationally cheap framework where the analyticity of the quark propagator with a non-positive spectral 
function can be exploited on the one hand and high energy QCD constraints can be imposed. Nonethe- 
less, according to a revealing interview with Rev. Matthews, a former Cannon of the Durham Cathedral, 
we should ask not only how it works, but why ! 
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